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Abstract 

I demonstrate how confinement in Coulomb-gauge QCD makes quark-quark 
states of the color anti-triplet (diquarks) move out of the physical spectrum. Mesons 
as color singlet quark- antiquark states, on the other hand have finite masses and for 
highly excited states in the meson spectra effective restoration of chiral symmetry 
can be observed. 

In Coulomb-gauge QCD the 00-component of the gluon propagator D^ u {x,t) has an 
instantaneous part Vc(x)S(t) and confinement means that — Vc(x) — ► oo for \x\ — ► oo. It 
was shown that Coulomb confinement is a necessary condition for confinement, i.e., that 
the gauge invariant quark- antiquark potential Vw{r) goes to infinity for r — > oo [lj. An 
almost linearly rising Coulomb potential has been suggested [2], which was also confirmed 
by results from the lattice [3]. In momentum space it becomes 

V c (\k\) = ^ (1) 

where ac is the Coulomb string tension. Based on previous works [I] we performed a study 
of the mechanism of Coulomb confinement in the Dyson-Schwinger-Bethe-Salpeter frame- 
work [5] in Rainbow-ladder approximation [6J. We took into account only the Coulomb 
potential ([[]) and neglected transverse gluons and noninstantaneous contributions to D 00 . 
In that way all integrals over ko can be performed analytically and one has to deal with 
three-dimensional integral equations only. (CQ) has also an unrealistic ultraviolet (UV) be- 
havior. However, it has the advantage that it produces no UV divergences, thus making 
renormalization not necessary. Due to all these approximations some physics is lost and 
the model is not expected to provide realistic quantitative results but some qualitative 
insight into the physics of confinement in QCD. Since the axial-vector Ward-Takahashi 
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identity is satisfied, chiral symmetry and its dynamical breaking are respected. In partic- 
ular the pion mass becomes zero in the chiral limit, i.e., for vanishing current quark mass. 
The potential ([I]) causes infrared (IR) divergences which are regulated by introducing an 
IR regulator /ini and replacing k 2 — > k 2 + /ij R (here and in the following k = \k\). The IR 
limit is then taken by means of /iir — > 0. The essential point is that the integral 

A 775 J, ■ ..2 5, = (2) 



2n 2 J *((p-q\T + tf R f 2 m 
diverges in the IR limit and one can write 

2^5 / d3 1 V c(\P~ ?!)/(?) = ^ jd'q S(p- q)f(q) + IR finite term. (3) 
For the gap equation of the quark propagator S(p) the ansatz 

S-'ip) = -i (loPo - 7 " pC(p) - B{p)) (4) 
leads to a coupled system of two integral equations 

B(p) = m+ ±ftf qVo {W-$)^ (5) 

c( P ) = i + ^/AW-?1)p'?4t- ( 6 ) 



where p = p/p, m is the current-quark mass, io(p) = y M 2 (p) + p 2 and M(q) = ^4 is 
called the quark mass function. Using (j3J) yields 



with 



B(p) = ° c M{ f\ + IR finite term = ^L^M + IR finite term, (7) 

2fi m LO{p) 2/x m LO{p) 

C(p) = ° C _) . + IR finite term = + tr finite tenn ( 8 ) 

2/i m u;(p) 2fi 1R uj{p) 

u(p) = J B 2 (p) + p 2 C 2 {p)) = + IR finite term. (9) 

2/ilR 

The functions B{p) and C{p) diverge like /Ijr 1 but the mass function is IR finite. The IR 
behavior for M(p) and o>(p) for m = is demonstrated in the upper plots of Fig. 1, while 
in the lower plots the same quantities for constant //ir = 10~ 3 yfoc but different current 
quark masses are shown. The mass function converges to a finite function. For large p 
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Figure 1: Mass function M(p) (left plots) and u(p) (right plots): In the upper plots the 
IR behavior for current quark mass m = and in the lower plots the results for different 
current quark masses at constant //ir = 1CT 3 ^fd~c is shown. 

it goes to the current quark mass and for small p it gets a dynamical mass which is of 
approximately the same absolute size for different current quark masses. u;(p)//ni indeed 
becomes a constant which is independent of the current quark mass. 
The Bethe-Salpeter equation (BSE) for a meson with mass M is 



d 4 q 



Vc(\p- q\) loS(q + M/2, q) X (q, M)S(q - M/2, g) 7o . (10) 



For the pseudoscalar meson (pion with M = m f ) the Bethe-Salpeter amplitude \{Pi m -n 
contains three (pseudoscalar, axial-vector and tensor) components: 



X(p, m*) = P P G»)75 + rn n P A (p)7o75 + m n P T (p)p ■ 77075- 



IF 
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The BSE is reduced to a coupled system of integral equations 



<o(p)h(p) = ^/AW-^Ul + ^A (12) 



,2 



-M-^W - %) + ^/^ ca? _ 9l) «dp 9h) (13) 



for the two functions 



g{p) 



hip) = (14) 
w(p) 



%) + 2^-P A {p) + 2^f T P r (p) 



^ 2 (p) - T L q (p) "(pY 



(15) 



The IR behavior of Eqs. (I12II13P follows by using Eq. (jUj) on the left and Eq. ([3$ on the 
right hand sides, respectively which yields 

(Jq (Jq 

hyp) + IR finite term = hyp) + IR finite term, (16) 

2/i IR 

(Jq (Jq 

-g[p) + IR finite term = gyp) + IR finite term. (17) 



2/^m 2/irR l 

Obviously the IR divergences cancel in both equations and all physical observables, in 
particular the mass, are determined by the IR finite terms. The functions hyp) and g(p) 
have an IR finite limit, too. This is demonstrated for m = in the left plot of Fig. 2. 

For two quarks in the SU(3)c anti-triplet state (diquark) a color factor | enters into 
the BSE kernel. Apart from that for a scalar diquark with mass msr> one has the same 
integral equations as for the pion but the IR divergent terms on the right hand sides 
of Eqs. (I16II17P are reduced to only half their sizes. Thus there must be additional IR 
divergences. Replacing (11 6 II 171) by 

2 

° C -hyp) + IR finite term = -p-h(p) + ^^-gyp) + IR finite term, (18) 



2/xir 4/xi R 

2 

° g(p) SD ^ m 9(p) + IR- finite term = hyp) + — gyp) + IR finite term (19) 



2/xir 2a c 4/xir 

and solving for msD and gyp) yields 



msD = tt^ ; g(p) = -^h(p). (20) 
2/im o-q 
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Figure 2: IR limit of the functions h and g for the pion (left) and the scalar diquark 
(right) for m = 0. In all cases the normalization has been chosen such that h(0) = 1. 



The divergences are now balanced by introducing a relation between h(p) and g(p) and 
making the diquark mass IR divergent. In that way the diquarks are removed from 
the physical spectrum. This mechanism of confinement applies not only for the scalar 
diquark but for diquarks of all quantum numbers. Numerically we have reproduced the 
IR divergence of the mass for scalar and axial- vector diquarks [6]. On the other hand, 
the shape of the functions h(p) and g(p) converges in the IR limit, which is demonstrated 
in the right plot of Fig. 2. For that reason not only mesons but also diquarks have IR 
finite radii. Calculations for the electromagnetic form factor and the charge radius of the 
the pion in Coulomb-gauge QCD have already been performed earlier [TJ. Explicitly, the 
charge radius of the pion is given by 
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u 2 



d?p f 3 
(2tt) 3 \ 32cj 4 (p) 



2uj 2 (p) + (m( P ) -2p 2 M' 



g{p)h(p) 



W(p)h(p) + g( P )h'(p)} + ^ [g(p)h"(p) + g"(p)h(p) - 2g'(p)ti(p)]\ + 0( m ) (21) 



with 

and of the scalar diquark by 
3/im f d 3 p 
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(2tt) 3 [ Qu 4 (p) 



d 3 p 



g(p)h(p), 



2lu 2 ( P ) + (M(p) - 2p 2 M\p) 



[22) 



h 2 {p) 
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+2h(p)h'(p) + ^ [h(p)h"(p) - h'\p)] | + OW (23) 

with 

A/? D = 48 m / 4^ 2 W- ( 24 ) 



(2*) 



In both cases ' means ^?y- For quark mass m = we have obtained the results r n = 

4.3 cr^ 1 ^ 2 and tsd = 6.0 o- c ^ 2 [6]. Notice that A/^ converges to a finite value while A/f D 
goes to zero like /^ir. That means that for the diquark only the shape of h(p) converges 
but its size diverges like ■ Due to (!20l) . the size of g(p) goes to zero like /z^ 2 on the 
other hand. 

Finally I present results for the highly excited meson spectra in the chiral limit [5]. 
There are certain phenomenological evidences that in highly excited hadrons the chiral 
(SU(2)l x SU(2)r) and U(1)a symmetries are approximately restored (for a review see 
[9]). The states fall into approximate multiplets of SU(2)l x SU(2)r and the mass 
splittings within the multiplets vanish at radial quantum number n — > oo and/or spin 
J — > oo. Furthermore the splittings within a multiplet become much smaller than between 
the two subsequent multiplets. The reason for this "effective" symmetry restoration is 
that excited hadrons gradually decouple from the quark condensates due to a diminishing 
importance of quantum fluctuations [10]. I restrict the discussion here to scalar and 
pseudoscalar mesons. Given the complete set of standard quantum numbers /, J PC , the 
multiplets of SU{2) L x SU{2) R for J = are [II] 

(1/2, l/2) a : 1,0-+ M 0,0 ++ and (1/2, l/2) b : 1, 0++ <— > 0, (T+. 

The BSE for a scalar meson with mass m ++ is reduced to the coupled system of integral 
equations 

HpMp) = ^ / A ,1) w + (M,) + , (25) 




Hp) + 7^Jd 3 qV c (\p-q\)p-qg(q). (26) 



For highly excited states the typical momenta of the quarks become large. For large 
momenta, however, the mass function M(p) becomes small. Setting M(p) = in the 
second equation for the pseudoscalar meson ([TBI and the first equation for the scalar 
meson ( |25l) gives just the second equation for the scalar meson ( |26l) and the first equation 
for the pseudoscalar meson (1T2T) . respectively. For large momenta with M(p) m the two 
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Figure 3: Angular (left) and radial (right) Regge trajectories for isovector mesons. Mesons 
of the chiral multiplet (1/2, l/2) a are indicated by circles, of (1/2, 1/2)6 by triangles, and 
of (0, 1) © (1, 0) by squares ( J ++ and J for even and odd J, respectively) and diamonds 
(J and J ++ for even and odd J, respectively). 

systems of coupled integral equations become approximately the same. This can explain 
why pseudoscalar and scalar mesons with large n become approximately degenerate. For 
states with J > similar arguments hold but there are additional states which fall in the 
multiplets (0, 0) and (0, 1) © (1, 0). The Bethe-Salpeter amplitudes for mesons with large 
J become strongly suppressed for small momenta and already states with n = become 
approximately degenerate. In our model the numerical results for the meson spectra up to 
n, J = 6 show a very fast restoration of both SU(2) L x SU{2) R and U(1)a symmetries with 
increasing J and essentially more slow restoration with increasing n. The excited states 
lie on approximately linear radial and angular Regge trajectories which is demonstrated 
in Fig. 3. In the limit n —>■ oo and/or J — > oo one observes an approximate degeneracy of 
all states within the representation [(0, 1/2) © (1/2, 0)] x [(0, 1/2) © (1/2, 0)] that combines 
all possible chiral representations for systems of two massless quarks [TT] . 
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